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Abstract. In this paper, we revisit generalized Chaplygin gas (GCG) model as a unified dark matter 
and dark energy model. The energy density of GCG model is given as pccG / PGCGQ ~ [Bs + (1 — 
Bs)a-3(^+"']i/(^+"\ where a and are two model parameters which will be constrained by type la 
supernova as standard candles, baryon acoustic oscillation as standard rulers and the seventh year full 
WMAP data points. In this paper, we will not separate GCG into dark matter and dark energy parts 
any more as adopted in the literatures. By using Markov Chain Monte Carlo method, we find the result: 
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1 Introduction 



Since 1998, the type la supernova (SNIa) observations [T] 
have shown that our universe has entered into a phase of 
accelerating expansion. During these years from that time, 
many additional observational results, including current 
Cosmic Microwave Background (CMB) anisotropy mea- 
surement [2], and the data of the Large Scale Structure 
(LSS) [3], also strongly support this suggestion. And these 
cosmic observations indicates that baryon matter com- 
ponent is about 4% for total energy density, and about 
96% energy density in universe is invisible. Considering 
the four-dimensional standard cosmology, this accelerated 
expansion for universe predict that dark energy (DE) as 
an exotic component with negative pressure is filled in 
universe. And it is shown that DE takes up about two- 
thirds of the total energy density from cosmic observa- 
tions. The remaining one third is dark matter (DM). In 
theory mounting DE models have already been constructed, 
for the reviews please see [l[5Jini[71|51[Sl[in| . But there exists 
another possibility that the invisible energy component is 
a unified dark fluid, i.e. a mixture of dark matter and dark 
energy. Actually, based on the Einstein's gravity theory 
and the observed geometry structure of our universe, one 
can define the so-called dark fluid 
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dark 
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obs 



(1) 



where G^i/ is the Einstein tensor and T°^^ is the ob- 
served energy-momentum tensor. This property is dubbed 
as dark degeneracy [TTlfT^ . 

In these unified dark fluid models, the Chaplygin gas 
(CG) and its generalized model have been widely studied 
for interpreting the accelerating universe P^fHllTCllTB] . 
The most interesting property for this scenario is that, two 
unknown dark sections in universe-dark energy and dark 
matter can be unified by using an exotic equation of state. 
The original Chaplygin gas model can be obtained from 
the string Nambu-Goto action in the light cone coordinate 
|17) . For generalized Chaplygin gas (GCG), it emerges as 
an effective fluid of a generalized d-brane in a (d -|- 1, 1) 
space time, and its action can be written as a generalized 
Born-Infeld form [T4l . Considering that the application of 
string theory in principle is in very high energy when the 
quantum effects is important in early universe [17], the 
quantum cosmological studies of the CG and the GCG 
has been well investigated in Ref. [T7| and [T5] . 

The GCG model is characterized by two model pa- 
rameters Bs and a. To constrain the model parameter 
space, the GCG model has been confronted by cosmic ob- 
servations, please see 1!) 20 '22 for examples. In Ref. |19) . 
the geometric information from SN la, the baryon acous- 
tic oscillation (BAG) and shift parameter R, la{z^) and 

from cosmic microwave background radiation (CMB) 
were used, where Bs = 0.73j:^:^^ and a = -OmtoH was 
obtained. In these papers, the CCG was decomposed into 
two parts: effective dark matter and dark energy. And the 
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effective dimensionless matter energy density was given 
as a„^ = i76 + (1 - - where and a 

are dimensionless model parameters, for the details please 
see Eq. ^ of this paper. In Ref. [50], we obtained the re- 



sults: -Be 



0.7475;°;°i^ and a 



-0.0256 



-1-0.1760 
-0.1326' 



where 



high redshift Gamma ray bursts data points were added 
as 'standard candles' and the decomposition was also em- 
ployed. In these papers, the perturbation of GCG was not 
considered. So the values of model parameter a < was 
included. In the following sections of this paper, we will 
see that the values of a > are mandatory for the stabil- 
ities of GCG perturbations. As is pointed by the authors 
in |16) . when the phantom region is forbidden, one has 
a unique decomposition of GCG. However, the possibil- 
ity of phantom dark energy has not been ruled out by 
current cosmic observations ^I]. Currently there is still a 
lack of physics principle to do any decomposition. In fact 
one can list many kinds of decompostions, for example 
pc = PcoO'~^ defined as cold dark matter and the remain- 
ing part pde = Pgcg — Pc defined as dark energy. The 
serious problem would be that there does not exist any 
decomposition at all, because it is a whole energy compo- 
nent. So any kind of decomposition would be non-proper. 
In Ref. [22 , the authors took GCG as a whole energy com- 
ponent, where the perturbations of this unified dark fiuid 
was also discussed by WMAP-5 data and matter (baryon) 
power spectrum. In that paper, a tighter constraint was 
obtained when the matter (baryon) power spectrum was 
included. Based on these points, in this paper, we will 
use the currently available observational data sets: the full 
information of 7-year WMAP CMB data, BAG and the 
Umon2 SNIa data to constrain the GCG model as the uni- 
fication of dark matter and dark energy. We will see that 
when the full information of CMB is included, a tighter 
constraint will be obtained. 

This paper is organized as follows. In section II, the 
GCG model as the unification of dark matter and dark 
energy is introduced briefly. Based on the observational 
data, we constrain the GCG model in section III. Section 
IV is the summary. 



2 Main equations in generalized Chaplygin 
gas model 

The equation of state (EoS) of GCG model is given in the 
form of 

PGCG = -A/pqcg (2) 

where A and a are model parameters. In general, for a 
spatially non-fiat FRW universe, the metric is 



-df + a^{t) 



1 



1 — kr'^ 



dr^+r^ide^ + 



sm 



(3) 

where k = 0, ±1 is the three-dimensional curvature and a 
is the scale factor. Using the energy conservation of GCG, 
one can rewrite its energy density as 



PGCG = Pgcgo Bs + {1 - Bs)a 



-3(1+q)] 1+' 



(4) 



where Bs = A/ p]^^qq and a are model parameters. Form 
Eq. (|4|), one can find that < i?s < 1 is demanded to 
keep the positivity of energy density. If a = in Eq. 
the standard ylCDM model is recovered. Taking GCG as 
a unified component, one has the Friedmann equation 



= Hi (1 - r?;, - a - ^k) Bs + {1- B 



-3(l+a) 1 + 



} 



(5) 



where B is the Hubble parameter with its current value 
Bq — 100/ikm s~"'^Mpc~^, and J7i {i — b,r,k) are dimen- 
sionless energy parameters of baryon, radiation and effec- 
tive curvature density respectively. In this paper, we only 
consider the spatially flat FRW universe. 

To study the effects on CMB anisotropic power spec- 
trum, the perturbation evolution equations for GCG would 
be studied. We treat GCG as a unified dark fluid which 
interacts with the remaining matter purely through grav- 
ity. With assumption of pure adiabatic contribution to the 
perturbations, the speed of sound for GCG is 



-aw, 



2 ^Sp 
"'^ 5p P 
where w is the EoS of GCG in the form of 



Bs + (1 - S,)a-3(i+") ■ 



(6) 



(7) 



From the above equation, one can find that the values of 
w are non-positive. So to keep the non-negativity of sound 
of speed, a > is required. 

In the synchronous gauge, using the conservation of 
energy-momentum tensor T^^.^ = 0, one has the perturba- 
tion equations of density contrast and velocity divergence 
for GCG 

Sgcg = -(1 + w)i9GCG + ^)- in{cl - w)5gcg (8) 

OccG = -n{l - 3cI)9gcg + T^k^SacG - k^aacm 

following the notation of Ma and Bertschinger j23, . For the 
perturbation theory in gauge ready formalism, please see 
[M] . The shear perturbation ctgcg = is assumed and the 
adiabatic initial conditions are adopted in our calculation. 

In this kind of unified dark fluid matter model, the 
averaging problem is involved for the non-linear perturba- 
tions. This problem comes from the fact {p) ^ p{{p)) as 
pointed out by authors [3S] . Actually one can check that 



ij>)^-{A/p^^)^-A/{pr=p{{p)), 



(10) 



in the case of a 0. But for the linear case {5 — Sp/{p) -C 
1), one has [25j 



(p) = = ~A{p)-"{{l - aS)) = ~A{py 



(11) 



So we must be careful when we consider the large structure 
formation in this kind of unified dark fluid model. We are 
not going to discuss this problem deeply, because the main 
task of this paper is to discuss its effects on CMB and to 
constrain the model parameter space by using SN la, BAG 
and CMB data points. 
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3 Constraint method and results 
3.1 Method and data points 

To constrain the parameter space, we use Markov Chain 
Monte Carlo (MCMC) method which is contained in a 
publicly available cosmoMC package including the 
CAMB ^27j code to calculate the theoretical CMB power 
spectrum. We modified the code for the GCG as a unified 
fluid model with its perturbations included. The following 
7-dimensional parameter space is adopted 



Prameters 



Mean with errors 



P = {cjb, 0s, T, a, Bs,n,,\og[W^As]} 



(12) 



where Wb = fl^h? is the physical baryon density, Os (mul- 
tiplied by 100) is the ration of the sound horizon and angu- 
lar diameter distance, t is the optical depth, a and are 
two newly added model parameters related to GCG, Us is 
scalar spectral index, Ag is the amplitude of of the initial 
power spectrum. Please notice that the current dimension- 
less energy density of GCG ficcG is a derived parameter 
in a spatially flat [k = 0) FRW universe. So, it is not in- 
cluded in the model parameter space P. The pivot scale 
of the initial scalar power spectrum kso — 0.05Mpc^^ is 
used. We take the following priors to model parameters: 
Wfc e [0.005, 0.1], 0s e [0.5, 10], r e [0.01, 0.8], a £ [0, 0.1], 
Bs e [0, 1], e [0.5, 1.5] and log[10i%] £ [2.7,4]. In ad- 
dition, the hard coded prior on the comic age lOGyr < 
to < 20Gyr is imposed. Also, the weak Gaussian prior 
on the physical baryon density oji, — 0.022 ± 0.002 [55] 
from big bang nucleosynthesis and new Hubble constant 
Ho = 74.2 ± 3.6kms"^Mpc"^ [29] are adopted. 

The total likelihood C oc is calculated to get 

the distribution, here is given as 



e 

T 

a 
Bs 

ris 

log[10i°A«] 

^GCG 

Age/Gyr 
Ho 



n n00Q+U-LlL)U548+U.UUiiU 
"•U^^y-0. 000545-0. 00104 
1 r,-og+0. 00248+0. 00483 
i.Uy:00_o. 00248-0. 00496 
n nccQ+0. 00656+0. 0249 
U.UOOy_o 00713-0.0241 
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n Qcn+0-0144+0.0311 
U.yOU_o 0145-0.0274 
o r|Q9Q+0. 0337+0. 0679 
O.UOZy_o 0341-0.0678 
n Qcc+0. 00166+0. 00319 
U.»00_o. 00164-0. 00336 
1 C77+0. 110+0. 218 
iO.O/ / -0.112-0.218 
n n/l/IR1+0 00164+0.00336 
U.U1'±U±_o 00166-0.00319 
in /lQS+1.200+2.321 
iU.4y8_i 207-2.454 
71 790+1-535+3.112 
' ^- ' ^^-1.525-3.0408 



Table 1. The mean values of model parameters with la and 
2(7 errors from the combination WMAP+BAO+SN. 



where Cs = 1 / ^3(1 
substitution [Ml 



R) is the sound speed. We use the 



d^— — Trs{zd), 



' fs{zd) 



(16) 



to obtain unbiased parameter and error estimates, where 
dz = rs{zd)/ Dv(z), fs is evaluated for the fiducial cos- 
mology of Ref. [31], and Zd is obtained by using the fit- 
ting formula [31] for the fiducial cosmology. Here Dv{z) = 
[(1 -I- z)'^D\cz/H{z)]^/^ is the 'volume distance' with the 
angular diameter distance Da- The 557 Union2 data with 
systematic errors are also included [31]. For the detailed 
description of SN, please see Refs. [55] . 



X 



XCMB 



Xbao 



XSN- 



(23) 3.2 Fitting Results and Discussion 



The CMB data include temperature and polarization power 
spectrum from WMAP 7-year data [50] as dynamic con- 
straint. The geometric constraint comes from standard 
ruler BAO and standard candle SN la. For BAO, the 
values {rs{zd)/DY{0.2),rs{zd)/Dv{0-5)} and their inverse 
covariant matrix [31] are used. To use the BAO informa- 
tion, one needs to know the sound horizon at the redshift 
of drag epoch Zd- Usually, Zd is obtained by using the 
accurate fitting formula [32] which is valid if the matter 
scalings pb oc and pc oc are respected. Obviously, 
it is not true in our case. So, we find Zd numerically from 
the following integration |55] 



rvo 

T{m)^ / di-j'Td 

^^dl] Xe{z)(JT 

da R 



(14) 



where R — iph/Ap^^ ut is the Thomson cross-section and 
Xe{z) is the fraction of free electrons. Then the sound hori- 
zon is 

rs{zd) = / d7]Cs{l + z). (15) 



We generate 8 independent chains in parallel and stop 
sampling by checking the worst e- values [the variance(mean) /mean("v 
of 1/2 chains] i? — 1 of the order 0.01. The calculated re- 
sults of the model parameters and derived parameters are 
shown in Table. [TJ where the mean values with la and 
2ct regions from the combination WMAP-I-BAO-I-SN are 
listed. The minimum is 8010.420 which is larger than 
that xlitn = 8009.116 for ylCDM model with the same 
data sets combination. Correspondingly, the contour plots 
are shown in Figure [TJ 



From the Table [T] and Figure [1] one can clearly see that 
a tight constraint is obtained when the full information of 
CMB data is included. For the small values of a, one finds 
that GCG model is very close to yiCDM model. 

To understand the effects of model parameters to the 
CMB anisotropic power spectra, we plot the Figure [U 
where one of two model parameters a and Bg varies, where 
the other relevant parameters are fixed to their mean val- 
ues as listed in Table [T] The panels of Figure [2] show the 
effect of parameter a and B to CMB power spectra re- 
spectively. The model parameters a modifies the power 
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Fig. 1. The ID marginalized distribution on individual param- 
eter and 2D contours with 68% C.L. and 95% C.L. by using 
CMB+BAO+SN data points. 

law of the energy density of GCG, then it makes the grav- 
ity potential evolution at late epoch of the universe. As 
results, one can see Integrated SachsWolfe (ISW) effect 
on the large scale as shown in the left panel of Figure [2J 
In the early epoch, GCG behaves like cold dark matter 
with almost zero EoS and speed of sound cj, therefore the 
variation of the values of a will change the ratio of energy 
densities of the effective cold dark matter and baryons. 
One can read the corresponding effects from the variation 
of the first and the second peaks of CMB power spectra. 
The parameter Bs is related to the dimensionless density 
parameter of effective cold dark matter J7c0- Decreasing 
the values of B^, which is equivalent to increase the value 
of effective dimensionless energy density of cold dark mat- 
ter, will make the equality of matter and radiation earlier, 
therefore the sound horizon is decreased. As a result, the 
first peak is depressed. The Figure [3] shows CMB power 
spectra with mean values listed in Table [T] for GCG and 
ylCDM model, where the black dots with error bars de- 
note the observed data with their corresponding uncer- 
tainties from WMAP 7-year results, the red solid line is 
for GCG with mean values as shown in Table [1] the blue 
dashed line is for ylCDM model with mean values for the 
same data points combination. And the green doted line is 
for ylCDM model with mean values taken from [21] with 
WMAP-|-BAO-|-iJo constraint results. One can see that 
GCG can match observational data points and ylCDM 
model well. 




Fig. 2. The CMB C; power spectrum v.s. multiple moment 
I. The panels show the effects of model parameters a and Bs 
to CMB temperature anisotropic power spectra respectively, 
in each case the other relevant model parameters are fixed to 
the mean values as listed in Table [T] 

and sound speed at early epoch (a < 0.2). The almost zero 
sound speed is important for large scale structure forma- 
tion. One can also read that it behaves like dark energy 
with EoS w < at late time, which pushes the universe 
into an accelerated phase. Here we do not consider the 
large structure formation in this unified GCG framework. 
In fact, it is not a easy issue. We leave it for the future 
work. 



4 Summary 



Now, we'd like to discuss the properties of GCG in the 
framework of unified dark fluid way, because the interpre- 
tation in the decomposition way is strongly decomposition 
dependent. The physical properties of GCG are described 
by its EoS w and adiabatic sound speed which are de- 
termined by model parameters Bs and a. We can show 
the evolutions of w and with respect to scale factor a 
in Figure U) From the left panel of Figure HI one see that 
GCG behaves like cold dark matter with almost zero EoS 



In summary, we perform a global fitting on GCG model, 
which is treated as a unified dark matter and dark energy 
model, by using MCMC method with the combination of 
the full CMB, BAO and SN la data points. As a contrast 
to the reports in the literatures, we take GCG as an entire 
energy component and without any decomposition. Tight 
constraint is obtained as shown in Table [T] and Figure 
[T] The GCG model can match observational data points 
and ylCDM model well. For the small values of model 
parameter a, one can conclude that MCG model is very 
close to ylCDM model. And currently available data sets 
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Fig. 3. The CMB Cf^ power spectrum v.s. multiple moment 
I, where the black dots with error bars denote the observed 
data with their corresponding uncertainties from WMAP 7- 
year results, the red solid line is for the unified dark fiuid model 
with mean values as shown in Table [TJ the blue dashed line is 
for ylCDM model with mean values for the same data points 
combination. And the green doted line is for ylCDM model 
with mean values taken from ^ with WMAP-|-BAO-|-//o con- 
straint results. 



of CMB, BAG and SN la can not distinguish GCG model 
from ylCDM model. 
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